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ABSTRACT 

The detection of radio pulsars within the central few parsecs of the Galaxy would provide a unique 
probe of the gravitational and magneto-ionic environments in the Galactic center (GC) and, if close 
enough to Sgr A*, precise tests of general relativity in the strong- field regime. While it is difficult to 
find pulsars at radio wavelengths because of interstellar scattering, the payoff from detailed timing 
of pulsars in the GC warrants a concerted effort. To motivate pulsar surveys and help define search 
parameters for them, we constrain the pulsar number and spatial distribution using a wide range of 
multiwavelength measurements. These include the five known radio pulsars within 15' of Sgr A*, non- 
detections in high-frequency pulsar surveys of the central parsec, radio and gamma-ray measurements 
of diffuse emission, a catalog of radio point sources from an imaging survey, infrared observations of 
massive star populations in the central few parsecs, candidate pulsar wind nebulae in the inner 20 pc, 
and estimates of the core-collapse supernova rate based on X-ray measurements. We find that under 
current observational constraints, the inner parsec of the Galaxy could harbor as many as ^10'^ active 
radio pulsars that are beamed toward Earth. Such a large population would distort the low-frequency 
measurements of both the intrinsic spectrum of Sgr A* and the free-free absorption along the line of 
sight of Sgr A*. 



Subject headings: Galaxy: center — pulsars: general 



1. INTRODUCTION 

The discovery of one or more pulsars in the inner par- 
secs around Sgr A*, the massive black hole (MBH) at 
the center of our Galaxy, would provide an invaluable 
tool for studying the innermost regions of the Galactic 
center (GC). Most of the current understanding of the 
inner parsec comes from infrared observati ons of the nu- 
clear star cluster (for a recent review, see iGenzel et al.l 
|201CI[) . The nuclear star cluster is centered on Sgr A* 
and consists of young massive stars at a projected radius 
of r « 0.5 pc and a dense collection of B-stars (the "S- 
stars") within r < 0.04 pc with the closest orbit passing 
just 6 X 10-^ PC (r^IOO AU) from Sgr A* (ISchodel et al.l 
120021 : iGhez et al.ll2003| ). Two decades of monitoring the 
orbits of these S-stars has yielded the mass of the central 
object to be M = 4 x lO^Mtn, unambiguously classify ing 
it as a MBH (IGhez et al.ll2008t IGillessen et aI1[200l . 

Despite the success of tracking stellar orbits in the in- 
frared, the sensitivity of this method is ultimately limited 
by source confusion. The detection of a radio pulsar at a 
similar distance with an orbital period of Porb ^ 100 yr 
would provide unparalleled tests of gravity in the strong- 
field regime. The timing of such a pulsar could allow the 
measur ement of the spin o r quadrupole moment of the 
MBH (IPfahl fc Loebl 200l ILaeima fc WolszczanI UWA 



iWex fc Kopeikinlll99alLiu et al.ll2012D . Additionallv. a 
pulsar found anywhere in the inner few parsecs of the 
Galaxy would provide a useful probe of the GC environ- 
ment. The mere detection of a pulsar would place con- 
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straints on the star formation history (SFH) and mea- 
surements of the dispersion measure and pulse broad- 
ening times would provide information on the electron 
density distribution of the region. 

However, even with th e detection of almost 2 000 radio 
pulsars in the Galaxy ([Manchester et ani2005() and sev- 
eral directed searches of the GC, only five pulsars have 
been found w ithin 15' of Sg r A* and the closest of these 
is 11' away (iDeneva et al.l l2009t iJohnston et"all 120061 : 
iBates et al.ll2011| ). While these few objects indicate the 
existence of a GC pulsar population, the perceived dearth 
of pulsars near Sgr A* is the result of interstellar scat- 
tering from turbulent plasma, which temporally broad- 
ens pulses to approximately 20001/^^^ s (where z^ghz is 
the obs erving frequency in G Hz) at the center of the 
Galaxy (jCordes fc Lazicill2002f ). Pulse broadening makes 
it almost impossible to detect even long-period pulsars 
in periodicity searches at commonly used frequencies 
(j/ ~ 1 GHz). To mitigate the deleterious effects of inter- 
stellar scattering, periodicity searches of the GC have 
migrated to higher frequencies (v ~ 10 GHz). How- 
ever, since pulsars have power-law spectra of the form 
S^v) oc 1'°' (with a < 0), increasing the observing fre- 
quency also decreases the observable flux density. To 
date, high-frequency searches have produced no new de- 
tectio ns using existing 100 m class telescopes (|Deneval 
IMoHMacquart et al.ll2010l ). 

Even though the absence of pulsar detections in the cen- 
tral parsecs of the GC is well explained by scattering ef- 
fects, the existence of a GC pu lsar population was estab- 
lished by IDeneva et al.l ()2009f ) based on the five pulsars 
on the outskirts of the region that cannot be explained 
as foreground disk objects. Since future surveys can ben- 
efit from better knowledge of the pulsar populations in 
the GC, we use a suite of multiwavelength observations 
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Figure 1. Views of tlie inner GC region at radio, X-ray and in- 
frared wavelengths. From top to bottom, the GC is shown in ra- 
dio at 8.5 GHz as observed with the Green Bank Telescope (data 
courtesy Casey Law), in 0.5-7 keV X-rays as observed with the 
Chandra A CIS-I instrument, and at J band (1.25 fim) as observed 
in 2MASS IjSkrutskie et al.ll2006I V In each panel, a circle with ra- 
dius 25" (corresponding to 1 pc at 8.5 kpc) is centered on Sgr A* 
wi th J2000 coordinate s of (1 7^45'°40!0409, -29°00'28'.' 118) given 
bv lReid fc Brunthalerl |(20M '). 



to set constraints on the number and distribution of pul- 
sars in the inner regions of the Galaxy on ~100 pc and 
~1 pc scales. An illustration of the structure of the GC 
on these scales is shown in Figure [TJ 
In this paper, we present observational constraints on the 
pulsar populations in the GC. A brief overview of the con- 
ventions and notations used in this paper are presented 
in Section [21 In Section [3l population limits are set from 
the detections of pulsars in the inner 15' and from non- 
detections in the vicinity of Sgr A* on parsec scales. In 
Section |4l a catalog of steep-spectrum radio sources in 
the inner 150 pc is considered. Interferometric measure- 



ments of the spectrum of Sgr A* are used in Section \5\ 
to set upper limits on the pulsar population in the GC 
on arcsecond scales. In Section [6l Fermi observations of 
the diffuse gamma-ray flux of the inner degree of the GC 
are used to estimate the millisecond pulsar population 
in the GC. In Section [71 infrared observations of young 
massive stars are used to estimate the number of neutron 
stars produced in the inner parsec of the Galaxy. Chan- 
dra X-ray observations of pulsar wind nebulae are used 
to constrain the pulsar population in the inner 20 pc in 
Section[8l In Section[9l limits are set on the intrinsic neu- 
tron star population in the GC based on the estimated 
supernova rate. Finally, in Section [TUl the estimates are 
summarized and discussed. 

2. CONVENTIONS AND NOTATIONS 

We note briefly a few conventions and notations that we 
adopt in this paper. The name "Sgr A*" is used to de- 
scribe both the MBH and the compact radio source at the 
dynamical center of the Galaxy. In the event where a dis- 
tinction must be made (for example, Section[5]), "Sgr A*" 
is taken to mean the observed radio source. All distances 
from Sgr A* are given as a projected distance unless ex- 
plicitly stated otherwise. The distan ce from Earth to 
Sgr A* is taken to be d = 8.5 kpc (jChez et all 120081 : 
Gillcs scn ct al„ 2009). 

Pulse broadening times are taken from the N E2001 
electron density model of iCordes fc Laziol (|2002( ). The 
NE2001 model accounts for the geometry of the scatter- 
ing region in the GC, and as such, it is prefera ble over em- 
pirica l fits to pulsars in the Galactic disk (e.g. iBhat et aLl 
l2004f ). which do not consider the pa rticular scatte ring ge- 
ometry of the GC. As discussed in lDeneva et al.l (pOOQi) . 
NE2001 tends to overestimate the scattering times for the 
five pulsars closest to Sgr A* by factors of ~10^ — 10"^. 
However, these pulsars likely lie along the edges of the 
GC scattering region where slight (^0.1 kpc) changes to 
the line of sight distance can cause dramatic (factors of 
^10^) changes in scattering times with only modest (fac- 
tor of ~2) changes in dispersion measure. Additionally, 
we note that the scattering times of pulsars deeper in the 
GC will be highly constrained by the measured angular 
scattering of Sgr A* itself. 

Pulse broadening times are found at arbitrary observ- 
ing frequencies by scaling the 1 GHz value s given by 
NE 2001 as oc t/~^ ([Lambert fc RickettI 119990 . However, 
see lLohmer et al] ([2001!) for potentially significant devi- 
ations from this scaling in highly scattered pulsars. 
All pulsar population estimates are given as the num- 
ber of active radio pulsars beamed toward the Earth, 
where a beaming fraction of /{, = 0.2 is assumed for all 
pulsars. The fixed beaming fraction of ft = 0.2 comes 
from a simple model in which the magnetic dipole mo- 
ment is oriented randomly with re spect to the rotation 
axis (jEmmering fc Chevalie^ 119890 . A better empiri- 
cal fit to the data is provided by the period-dependent 
model of iTauris fc Manchested (|19980 , which typically 
finds fh ^ 0.1 for isolated pulsars. However, since the 
iTauris fc Manchested (|1998| ) model does not include mil- 
lisecond pulsars tha t may have beaming fractions as high 
as /{, « 0.5 - 0.9 ([Kramer et al.lll9980 . we adopt the 
fb = 0.2 value as a population-wide representative value. 
As this factor only shows up in our analysis as a mul- 
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tiplicative constant, it is trivial to scale our results to 
different beaming fractions. 

In many of the limits presented below, it will be necessary 
to utilize a distribution for pulsar pseudo-luminosities 
{L = SdP). We ad opt the' power-law distribution of 
iLorimer et aD ()2006f ) with a minimum cutoff in prefer- 
ence to distributions t hat do not require cutoffs l ike th e 
log-normal model of iFaucher- Giguere fc KasDil (l2006h . 
Both arc two parameter models and the power-law distri- 
bution provides a much bett er empirical fit to the obser- 
vations (see, e.g.. Figure 6 of ILorimer et al.|[2006[ ). How- 
ever, if the log-normal distribution were used instead of 
the power-law in the calculations below, the total pul- 
sar populations predicted would be larger by factors of 
10-100 owing to the much larger fraction of low luminos- 
ity objects. As a result, our adoption of the power-law 
pseudo-luminosity distribution is a conservative one in 
the sense that using another distribution would predict 
a larger pulsar population. 

Lastly, we note that in some cases our constraints will 
involve all types of pulsars, while others involve one of 
two subsets of pulsars: the "canonical" pulsars (CPs) 
and recycled or "millisecond" pulsars (MSPs) . CPs have 
periods of P ~ 1 s, surface magnetic fields B ~ 10^^ G 
and active radio lifetimes of r ~ lO'' years, while MSPs 
have periods -P ^ 10 ms, low surface magnetic fields B < 
10^ G and active radio lifetimes of t '-^ 10^ — 10^" years. 

3. CONSTRAINTS FROM PULSAR SURVEYS OF THE GC 

Several pulsar searches have been conducted in the inner 
degree of the Galaxy. To date, only five pulsars have been 
detected within 15' of Sg r A*, with none closer than 11' 
(jManchester et al]|2005lf . To have any chance of mak- 
ing detections in the inner few arcminutes of the Galaxy, 
higher observing frequencies must be used to overcome 
the roughly 2000i^Qj![2 s broadening times caused by scat- 
tering. Deep searches have been attempted at frequencies 
from 4 to 15 GHz, but have mad e no detections in the 
inner few pa.rsecs around Sgr A* (iJohnston et all [20061: 
lDeneval[2010l: iMacquart et al.ll2010D . Both the pulsar de- 
tections in "low-frequency" (2-3 GHz) surveys and the 
absence of detections in high-frequency (4-15 GHz) di- 
rected searches can be used to constrain the GC pulsar 
population. 

3.1. GC Pulsar Detections in Low- Frequency Surveys 

The five known pulsars within 15' of Sgr A* currently 
provide the best direct evidence for an intrinsic GC 
pulsar population. Of these five pulsars, two were 
detected at 3.1 GHz with the Parkes radio telescope 
([Johnston et al.l[2006l ) and three^ were detected at 2 GHz 
with the Green Bank Telescope ([Deneva eFall [20091) . In 
each of these surveys, the expected number of detectable 
disk pulsars in the field of view is <C 1. Thus, the de- 
tections strongly suggest a pulsar population in the GC 
that is distinct from that of the disk. 
In an attempt to constrai n the number and sp atial distri- 
bution of the GC pulsars, [beneva et al.l ([2009| ) simulated 
the pulsar population to determine what would be consis- 
tent with the survey detections. A simple two-component 

^ On e of these pulsars, J1746— 2850, was independently discov- 
ered bv lBates etahl l l20Tll) in a Parkes 6.5 GHz multibeam survey. 



density model of the form 

nGcocexp(^--A_^exp(^--^^ (1) 

was adopted with Hqc = 26 pc fixed to coincide with 
the scale height of the sc attering screen in the N E2001 
electron density model of iCordes fc Laziol (|2002D . The 
density distribution was normalized so that, for a given 
i?(3C; there were a total of Nqq pulsars (not necessarily 
beamed toward Earth) associated with the population. 

Using Monte Carlo methods, iDeneva et al.l ([20091 ) then 
generated 1000 pulsar populations consistent with Equa- 
tion [T] for each {Rqc,Nqc) pair and determined how 
many of these pulsars would have been detected in their 
2.1 GHz survey. Searching over a grid of values and 
performing a maximum likelihood analvsis. IDeneva et al.l 
pOOa) found lower bounds of Ngc > 2000 and Rgc > 
0.3 kpc for the parameters. Thus, this analysis provides 
additional evidence for an intrinsic pulsar population in 
the GC. 

Although IDeneva et al.l ([2009f ) have set a lower bound on 
Ngci this parameter applies to the entire GC population 
and does not necessarily translate into a lower bound 
on the number of pulsars within a particular distance 
from Sgr A*. For example, a distribution with Nqc = 
2000 and Rgc = 0.3 kpc will produce a very different 
number of pulsars within 100 pc of Sgr A* than will 
a distribution with Ngc — 2000 and Rcc ~ 0.6 kpc. 
Regardless, the existence of a pulsar population in the 
GC is firmly established. 

We no te briefly that the primary reason the lDeneva et al.l 
([2009D analysis does not find upper bounds on the pa- 
rameters Ngc a-iid Rcc even after extending the grid 
to Ngc = 10^ and Rcc = 5 kpc is that the only con- 
straints come from the detections in the survey region 
(r < 50 pc). Incorporating survey results from the in- 
ner few degrees of the Galaxy would certainly introduce 
upper bounds to the parameters. 

With this in mind, it is instructive to consider the results 
at a fixed Rcc- From Figure 4 of IDeneva et"!!! (|2009j ). 
wc see that a wide range of Ngc values with Ngc ^ 
500 are equally likely for Rcc ~ 0.1 kpc. Thus, a very 
conservative lower bound on the number of pulsars in 
the inner 100 pc of the Galaxy that are beamed toward 
Earth is Npsr ^100 (where we have assumed a beaming 
fraction of 0.2). 

3.2. High- Frequency Pulsar Searches of the Central 
Parsec 

Recently, searches of the central few parsecs around 
Sgr A* have been conducted with t he Green Ban k Tele- 
scope ( GBT) at 5 GHz and 9 GHz ([Denevall2010[ ) and at 
15 GHz ([Macauart et al.l2010l ). No pulsars were detected 
in any of these searches . We follow a similar analysis to 
that of IMacquart et all ([2010) to estimate an upper limit 
to the pulsar population based on the absence of detec- 
tions. 

3.2. L Observations 

The 5 and 9 GHz observations were carried out by 
IDeneval ((Moh in 2006. Since no pulsar candidates were 
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detected, limits on the flux density of periodic signals 
may be set using the radiometer equation 



(2) 



where T^ys is the system temperature, G is the telescope 
gain, r] w 0.8 is a correction factor that accounts for 
system imperfections and the digitization of the signal, 
Nh = 16 is the maximum number of harmonics summed 
in a periodicity search, Npoi = 2 is the number of po- 
larization channels summed, Ai/ = 800 MHz is the re- 
ceiver bandwidth and Tots = 6.5 hr is the observation 
time. The telescope gain is given by G = 1.85 K Jy~^ 
and G = 1.8 K Jy~^ for observing frequencies of 5 and 
9 GHz, respectively. The value of m is determined by the 
detection significance threshold set at ma. In this FFT 
search the threshold was set to 6cr, so m = 6. 
The system temperature of the telescope is given by^ 
Tgys = Trec + Tbg. The receiver temperature of the GBT^ 
is 18 K at 5 GHz and 27 K at 9 GHz. The dominant con- 
tribution to the background temperature is the bright 
extended Sgr A Complex (comprised of Sgr A East and 
Sgr A West), which surrounds Sgr A*. We may set lower 
bounds on this b ackground using data from a multiwave- 
length survey bv lLaw et al.l (|2008f ). which imaged the GC 
at 1.4, 5 and 9 GHz. Only lower bounds may be set since 
a non-trivial iterative scheme was used to su btract out 
the no ise contributions from the atmosphere. iLaw et ahl 
(|2008| ) found the flux density of the Sgr A Complex 
to be 85 Jy beam~^ at 5 GHz and 39 Jy beam^^ at 
9 GHz, which translate to background temperatures of 
Tbg = 157 K and Tt,g = 70 K, respectively. Thus, the 
system temperature is Tsys = 175 K at 5 GHz and 
Tsys = 97 K at 9 GHz. From Equation [2l we find that 
Smin,v = 29 /xJy at 5 GHz and Smin,v = 17 /iJy at 9 GHz. 
Likewise, iMacauart et al.l ()2010[ ) conducted a search for 
pulsars at 14. 8 and 14.4 GHz using the GET in 2006 
and 2008. The lMacquart et al.l (|2010l ) search used a IOct 
detection threshold and combined observations to get an 
effective observation time of Tabs ~ 9.75 hr. The system 
temperature was determined to be T^ys « 35 K by firing 
a noise diode on a calibrator source. No pulsars were 
detected and a lOa detection threshold fiux density of 
Smin,v ~ 10 /iJy is sct from the 14.4 GHz measurements. 

3.2.2. Upper Limits on Observable Pulsar Population 

Assuming that a given pulsar in the GC will be detected 
with some probability , binomial statistics can be used 
to find the maximum number of pulsars consistent with 
zero detections. The simplest way to determine pd is to 
set it equal to the fraction of pulsars bright enough to be 
seen in each survey when placed at Sgr A*. This fraction 
can be estimated from the 1.4 GHz pseudo- luminosity 
function, which is given by dN/d\ogL cx L^^ taken over 
a range of pseudo- luminosities from L„„;„ = 0.1 mJy kpc^ 
to Lmax ~ 10"* mJy kpc^. We adopt a value of /3 = —0.7 

^ We take Tree to include all non-astronomical contributions to 
the system temperature from the receiver, spillover effects, and the 
atmosphere. 

^ As provided in the GBT Proposer's Guide: 
http : //www . gb . nrao . edu/gbtprops/man/GBTpg . pdf 



as an average of the two fits found by iLorimer et al.l 
(l2006l) of over 1000 pulsars observed in the Parkes Multi- 
beam Survey. The detection probability is then given by 
Pd = /L(i > Ldet), where 



L 



det 



Sn 



1.4 GHz 



(3) 



and Smin,v is the minimum detectable flux density of 
the search. Using the Smin v values found by each search 
(see Table [U we find that Vd = (0.027, 0.020, 0.015) at 
(5,9,15) GHz. 

Given the above detection probabilities and the lack of 
any detections in the surveys, the upper limits to the 
number of pulsars (at 99% confidence level) are found to 
be iV = (170, 230, 299) for the 5, 9, and 15 GHz observa- 
tions, respectively. However, since a pulsar with a spin 
period less than the pulse broadening time would have 
a greatly reduced chance of being detected, the calcu- 
lated upper limits are for pulsars with P > Tsc- Using 
the 1 GHz scattering time from the NE2001 model of 
iCordes fc Laziol (|2002f l and scaUng (cx i/^**) to the appro- 
priate frequency, the scatter broadening times are found 
to be = (4.2, 0.44, 0.05) s at observing frequencies of 
V = (5,9,15) GHz. 

3.2.3. Implications for Total GC Pulsar Population 



The upper limits of Section 13.2.21 are only valid for pul- 
sars within certain period ranges. To make an estimate of 
the total number of pulsars, knowledge of the underlying 
pulsar period distribution is needed. Though the period 
distribution of pulsars in the inner parsecs of the GC 
is entirely unknown, a reasonable approximation would 
be to assume the same distribution as the local pulsar 
population (to reduce observational biases). From the 
ATNF pulsar catalog® (jManchester et al.l [20051) . we see 
that there are 88 pulsars within 1 kpc of Earth. Of 
these, 5 have P > 4.2 s, 38 have P > 0.44 s and 55 
have P > 50 ms. This gives /p(P > 4.2 s) = 0.06, 
fp{P > 0.44 s) = 0.43 and /p(P > 50 ms) = 0.63 for 
the fractions of pulsars with periods greater than the 
scatter broadening times. 

Assuming these values are representative of the GC pop- 
ulation, we can estimate upper bounds on the total 
number of pulsars (regardless of period) to be N,nax = 
N/fp. Applying these corrections to the estimates of 
Section [3221 

gives N„iax = (2830, 535, 475) for observing 
frequencies of = (5, 9, 15) GHz. 

The half-power beam width of the GBT is 6{i') w 
150"(i^/5 GHz)-i. At Sgr A*, the projected radu of 
these beams are r w (3.0, 1.5, 1.0) pc at 5, 9 and 15 GHz, 
respectively. Thus, we estimate that there are as many 
as Nmax < 2830 pulsars beamed toward Earth within 
3.0 pc of Sgr A*. Nmax < 535 within r sa 1.5 pc and 



N„ 



< 475 within r w 1 pc. 



The results of this section are sum marized i n Ta - 
ble [1] In addition to the surveys of IDeneval ([2010l ) 
and lMacQuart "eFall ([20T0l) . we also include for refer- 
ence the less sensitiv e G C pointings from surveys by 
[Johnston etldl ([20060 and [Bates etall ([20 lit ). 

^ http: //www. atnf . csiro . au/research/pulsar/psrcat/ 
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Table 1 

Pulsar Upper Limits of High-Frcqucncy Surveys of the Central Parsees 



V 






fUL > Laet) 


NiP > Tsc) 


fp(P > Tsc) 


N/fp 


r 


Ref 


(GHz) 


(MJy) 












(pc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


4.85 


29 


4.2 


0.027 


< 170 


0.06 


< 2830 


3.0 


(1) 


6.6 


550t 


1.28 


0.002 


< 2220 


0.17 


< 13000 


4.0 


(2) 


8.4 


200 


0.46 


0.003 


< 1380 


0.43 


< 3200 


3.0 


(3) 


8.50 


17 


0.44 


0.020 


< 230 


0.43 


< 535 


1.5 


(1) 


14.4 


10 


0.05 


0.015 


< 299 


0.63 


< 475 


1.0 


(4) 



Notes. Columns 1 and 2 give the center frequency and minimum detection threshold, 
respectively, of a given survey. Column 3 gives the scattering time at the survey cen- 
ter frequency of a pul se originating from the GC as calculated with the NE2001 model 
llCordes fc Lazij 120021) . Column 4 gives the fraction of pulsars bright enough to be de- 
tected at the GC with the survey sensitivity (see Equation [S]! . Column 5 gives 99% upper 
limit values to the number of pulsars that could potentially have been seen by the survey 
that is, have P > Tsc) - Column 6 gives the fraction of pulsars in the ATNF catalog 
Manchester et aI1l2005l) that have periods longer than the scattering time and Column 7 
gives an estimated upper bound for the total number of pulsars in the GC. Column 8 gives 
the projected radial distance of the bcamsizc at the distance of Sgr A* and Column 9 gives 
the reference from which the survey was taken. 

t The value provided for Smin here is 9 times that given in Section 2.2 of I Bates et al.l 

11201 II) , which we b elieve to be a c alcul ation error. 

References. (1) FDeneval [20TOl: (2) IBates et afl |20TII : (3) IJohnston eTall 1200^ : (4) 
liVlacauart et al.ll2010l 



3.2.4. Caveats 

A number of assumptions are made in the upper limit 
estimates of the previous sections, so it is important to 
consider what happens if the assumptions fail. The first 
assumption is that the probability of detecting a pulsar 
(with period large enough not to be smeared out by in- 
terstellar scattering) is equal to the fraction of pulsars 
with 1.4 GHz pseudo-luminosities large enough to be de- 
tected at the GC. Since the detection probability only 
considers the best-case sensitivity of the telescope and 
ignores any effects of radio-frequency interference at the 
telescope end or intermittency at the pulsar end, it is 
likely to be an overestimate. An overestimate of the de- 
tection probability would result in an underestimate of 
the upper bound on the number of pulsars. 
Another assumption is that the 1.4 GHz pseudo- 
luminosity distribution of pulsars in the GC is the same 
as those in the Galactic field. However, if the GC region 
contains a larger number of low-luminosity MSPs, then 
the detection probability will again be overestimated and 
the upper limit underestimated. One could also imag- 
ine the case where enough bright young pulsars exist in 
the GC as a result of recent star formation to skew the 
pseudo-luminosity distribution to higher luminosities. In 
that case, the detection probability would be an under- 
estimate and the upper limits an overestimate. 
We have also assumed that the upper limits are for pul- 
sars with P > Tsc- However, as the scattering time gets 
to be a significant fraction of the pulsar period, it will 
start to smear out the signal and reduce the number of 
detectable harmonics. Since this is a gradual process, 
it will likely make some pulsars undetectable even with 
P > Tsc ■ As a result, we would be overestimating the 
minimal detectable period, which would cause the upper 
bound to be an underestimate. 

Finally, we have assumed that the period distribution 
of the pulsars within 1 kpc of Earth is representative of 
the GC population. This certainly does not have to be 



the the star formation histories of the Galac- 

tic field and GC are likely to be different. For example, 
the GC could potentially have a much higher concen- 
tration of young pulsars and old MSPs as compared to 
the Galactic disk. The increased stellar encounter rate in 
the GC could favor MSP production and any recent star- 
burst would favor young pulsars. In both cases, the pe- 
riods would be biased low. Thus, the period distribution 
would be skewed lower than assumed and the fraction of 
pulsars with periods greater than a certain value will be 
overestimated. This will result in the upper bound being 
underestimated. 

Since most of the assumptions made tend to decrease the 
upper limits, our estimates are best interpreted as the 
most restrictive upper bounds to the pulsar population 
in the inner few parsees of the GC. 

4. RADIO POINT SOURCES 

M otivated by t he study of GC pulsar se arch methods 
by iCordes fc Laz io (1993), iLazio fc Corde s (2008) per- 
formed a VLA survey of compact radio sources in the 
inner degree of the Galaxy. Though pulsars cannot be 
identified by their pulsed emission in an imaging sur- 
vey, promising pulsar candidates may be found by look- 
ing for steep-spectrum sources with angular diameters 
consistent with the angular broadening of point sources 
caused by scattering at locations near Sgr A* (wl" at 
1 GHz). Of th e 170 compact radio sources cataloged, 
ILazio fc CordesI ()2008D estimate that the number of pul- 
sars included is of order ~10. Based on this survey, up- 
per limits to the pulsar population within 1° (^150 pc) 
of Sgr A* may be estimated. 

4.1. Observations 

The survey was conducted at observing frequencies of 
1.4 and 5 GHz with the VLA in the A configuration. A 
total of 13 fields arranged in a hexagonal grid covered the 
region of the GC out to roughly 1° (150 pc) from Sgr A* 
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(the half-power radius of the VLA primary beam is 15' 
at 1.4 GHz). The typical resolution for the survey was a 
synthesized beam size of 2'.'4 x 1'.'3. 
So urces were ident i fied u sing a method similar to that 
of iLazio k, CordesI (119981 ) . Essentially, a histogram of 
intensities was constructed from the image of the primary 
beam for each field. If the field just contained noise, the 
intensity histogram would be a Gaussian with a mean of 
zero and a standard deviation equal to the thermal noise 
of 0.05 mjy per synthesized beam. Sources could then 
be determined by looking for deviations from this noise- 
only histogram. In practice, the histogram was found 
to have larger tails than a Gaussian, with zero mean 
and a standard deviation of «0.5 mJy beam~^. Since 
the resolution is comparable to the scattering size of a 
point source at the distance of Sgr A*, a lOcr detection 
threshold of Sdet ~ 5 mJy was adopted for the survey. 

4.2. Pulsar Population Estimate 

Given that Nobs ~ 10 pulsars were likely observed, the 
total pulsar population in the survey region can be es- 
timated as Npsr ~ Nobs/ fh, where is the fraction of 
pulsars luminous enough to be detected at the distance 
of Sgr A*. Taking the survey detection threshold to be 
Sdet ~ 5 mJy at 1.4 GHz, a pulsar must have a 1.4 GHz 
pseudo-luminosity of at least L^et = 360 mJy kpc^ 
to be detected at the distance of Sgr A*. The frac- 
tion of pulsars with L > Ldet can be determined from 
the 1.4 GHz pulsar luminosity function, which has the 
form dN/d\ogL oc L-^. The range of pulsar pseudo- 
luminosities is taken from Lmm = 0.1 mJy kpc to 
Lmax = 10^ mJy kpc^ and the exponent in the distri- 
butio n function is taken to be (3 ~ 0.7 (jLorimer et al.l 
|2006() . From this distribution, the fraction of pulsars lu- 
minous enough to be detected is /i, = 3 x lO^'^. For 
Nobs ^10 pulsars detected in the survey, we expect a to- 
tal population of Npsr ^ 3000 pulsars within 1° (150 pc) 
of Sgr A*. 

Though Npsr ~ 3000 is the nominal population estimate 
from the survey, a broader range result s if the assump- 
tions do not exactly hold. For instance, ILazio fc CordesI 
(j2008D estimate that Nobs ^ 10 of the unidentified steep- 
spectrum point sources will ultimately turn out to be 
radio pulsars. However, this number could range from 
zero to about 30. If one takes Nobs = 30, repeating the 
above analysis gives a pulsar population of Npsr ^ lO**. 
Additionally, one may consider the case in which no pul- 
sars were detected. Despite a follow-up o bservatio i i of 15 
of the pulsar candidates in this survey bv lDeneval ()2010l ) 
with the GET, none of the candidates have to date been 
confirmed. Assuming zero detections, an analysis sim- 
ilar to that in Section 13.2! gives an upper limit to the 
pulsar population of Npsr < 1500 at a 99% confidence 
level. These two extremes illustrate that although the 
survey allows for an estimate of the pulsar population in 
the inner degree of about 3000, the actual number could 
be below 1500 or as high as 10**. As a result, we take 
Npsr ^ 10^ as a conservative upper bound. 
Finally, we note that although the survey covers the re- 
gion within 1° (150 pc) of Sgr A*, there will be reduced 
sensitivity in the field centered on Sgr A*. The reduced 
sensitivity is the result of increased background temper- 



atures and greater sidelobes from the extended structure 
of the inner GC. In addition, since the scatter-broadening 
of point sources in the vicinity of Sgr A* (wl") is com- 
parable to the resolution of the survey, source confusion 
may become important in the innermost regions of the 
GC. Therefore, this survey would be largely insensitive 
to a fairly compact population of pulsars in the inner 
tens of arcseconds around Sgr A*. 

5. RADIO SPECTRUM OF SGR A* 

In this section, limits are placed on the maximum allow- 
able number of pulsars in the inner parsecs of the GC 
based on radio interferometer observations of Sgr A* on 
arcsecond scales (1" « 0.04 pc at 8.5 kpc). Due to the 
finite resolution of interferometers and the broadening of 
angular diameters as a result of the interstellar scatter- 
ing of radio waves, the Sgr A* radio source is actually 
extended (wl" at 1 GHz). Flux measurements of Sgr A* 
will therefore include a contribution from a collection of 
pulsars, if such a population exists. Although these pul- 
sars will be unresolved, upper limits on the total popula- 
tion may be set based on the total fiux density of Sgr A* 
in a manner analogous to s imilar constraints plac ed on 
pulsars in globular clusters (iFruchter &: Gosj|1990f ). 
We consider a model in which the observed flux density 
of Sgr A* is actually the combination of two components. 
The first component is that due to radio emission from 
the immediate environment of the MBH itself, which we 
assume is described accurately by high frequency obser- 
vations where the pulsar component is negligible. The 
second component is that due to the population of pul- 
sars near the MBH. This pulsar component becomes im- 
portant at lower frequencies both because radio pulsars 
typically have steep spectra and the angular resolution of 
radio telescopes scales with frequency such that a larger 
region around the MBH is sampled at lower frequencies. 
By requiring that this model flux be consistent with ex- 
isting observations, constraints may be set on the max- 
imum number of pulsars allowed in the inner parsec of 
the Galaxy. 

5.1. Observations 

Two different measurements of the spectrum of Sgr A* 
over a wide range of frequ en cies are co n sidere d (|An et al.l 
l2005t iFalckeet all [1991 . lAn et al.l (l2005l) conducted 
simultaneous measurements of Sgr A* from 300 MHz 
to 43 GHz using the VL A (A-configuration) and the 
GMRT. lFalcke et al.l (|1998D made simuhaneous measure- 
ments of Sgr A* using the VLA (A-configuration), the 
Berkeley-IUinois-Maryland Array, the Nobeyama 45 m 
telescope, and the Institut de Radioastronomie Mil- 
limetrique (IRAM) 30 m telescope from 1.4 GHz to 
235 GHz. Both groups observed a broken power-law 
spectrum with a break around 10 GHz. Since the power- 
law spectrum of pulsars decreases with increasing fre- 
quency, a population of pulsars contributes significantly 
only at lower frequencies. As a result, we consider the 
spectrum of Sgr A* only below the break frequency at 
10 GHz. 

5.2. Spectral Model 

We model the measured flux density of the compact ra- 
dio source Sgr A* as the sum of contributions from a 
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collection of pulsars and a point source associated with 
the MBH attenuated by free-free absorption according 



The effective resolution of the VLA observations is then 



to 



Ssg,{v)-- 



X exp(-!^y/i/^). 



^bh 



apsr' 



(4) 



The emission from the immediate vicinity of the MBH, 
5'bh,i/o, is taken to be a point source with a power-law 
spectrum with spectral index abh- is the number of 

pulsars contained in the solid angle of the effective point- 
spread function or beam size, which depends strongly on 
frequency (see Section 15. 3p . We use a simplified scaling 
for the free-free absorption that ignores the frequency 
dependence of the Gaunt factor. The free-free absorp - 
tion factor has a turnover frequency Vf ()An et al.ll2005[ ). 
The flux density per pulsar, 5'psr,i/o, is taken to be the 
mean of the 1.4 GHz pulsar pseudo -luminosity distribu- 
tion given by dN/dlogL cx (jLorimer "etall[2006h 

with a lower cutoff of L„iin = 0.1 mJy kpc and an up- 
per cutoff of L,nax = 10^ mJy kpc^. The mean observed 
flux density at 1.4 GHz is 5psr,i.4 = 99 fJ.3y. Since pulsar 
radio flux scales with frequency as a simple power-law, 
we can scale our fiux density as Si, oc z^"?^''. We fix 
ttpsr = —1.7 as a nominal value for the pulsa r spectral 
index (jMaron et all 120001: iLorimer et all 119951 ). but also 
consider values of —1.0 and —2.5 to test any major spec- 
tral index dependence. 



5.3. Effective Angular Resolution 

The number of pulsars included in a flux measurement 
of Sgr A* can be written as an integral of the number of 
pulsars per unit solid angle. 



(5) 



The integral is over the effective solid angle flcsiiy) 
(7r/4) 6*^^(1/), where 



;6) 



is the effective resolution with the subscript "0" repre- 
senting values at 1 GHz and the frequencies are in GHz 
units. The effective resolution is the quadrature sum of 
the resolution of the interferometer and the angular ex- 
tent of Sgr A* caused by scattering. The scaling for the 
synthesized array beam (oc assumes a fixed array 
configuration and the scattering diameter Ogc (oc i' '^) 
scales in conformanc e to me asurements of Sgr A* and 
OH/IR masers fe.g.. IfVail et al.. 1994' ). Our treatment 
assumes that scattered images are circular whereas in 
fact some are elliptical, but given that we are making 
order of magnitude estimates of pulsar numbers, the dif- 
ferences are not important. 

The angular diameter of Sgr A* is dominated by inter- 
stellar scattering at low fr equencies, with an observed 
major axis of Osco = 1''2 (jBower et al.ll2006D . For the 
VLA in the A configuration and ignoring any effects of 
foreshortening, the half-power beam width is ^b,, = 1'.'95 
([Bridle, 198^. 



cfi(!^) = l'.'2iy-^ (1 + 2.6 v 



(7) 



The two contributions are equal at v = 0.6 GHz, so at 
frequencies lower than this the resolution is completely 
scattering dominated and the resolution solid angle scales 
steeply with frequency as 

Additionally, a single data point measured with the 
GMRT will be considered in our analysis, so a similar 
effectiv e resolutioii must be constructed for this tele- 
scope. IRov fc Raol ([200l measure the resolution to be 
11'.'4 X 7'.'6 at 620 MHz. Converting this ellipse to a cir- 
cle of equal area and scaling to 1 GHz gives — 5'.' 77 
for the GMRT. Combining this with scattering as above, 
gives 



fcff 



{i^) = V'.lv'"^ (1 + 23.1 v^) 



1/2 



(8) 



For the GMRT, the two components of the effective res- 
olution are equal a.t v = 0.2 GHz. 

5.4. Candidate Pulsar Distributions 

We choose three physically-motivated distributions as 
model pulsar populations. The distributions are illus- 
trated below in Figure [2] 

The first (Model A) assumes a constant number of pul- 
sars per unit solid angle, dup/dfl = constant, so the num- 
ber of pulsars scales as 



(9) 



where fti is the solid angle enclosing A^i pulsars. In 
the fitting below, ili is set so that iVi gives the num- 
ber of pulsars in the inner parsec. Referring to Equa- 
tion 131 it may be seen that when scattering dominates 
the effective resolution, the contribution to the unab- 
sorbed spectrum from pulsars increases very rapidly as 
iV^(i^)z^"p=' oc i'^^-'^. Free-free absorption attenuates 
much of the flux, thus allowing a significant pulsar pop- 
ulation to remain hidden in spectral measurements. 
In Model B, we assume that dup/dil oc $7"° '', cor- 
responding to the surface density scaling observed for 
Wolf-Rayet and 0- star populations in the inner parsec 
(jCenzel et al.llMol) . This yields 



(10) 



However, the observed populations o f Wolf-Rayet and O - 
stars have an inner cutoff at 6* « 1" (|Bartkoet al.l[20T0l ). 
This core may affect the pulsar population in many ways, 
but we shall just consider two here. In both distributions, 
the pulsar surface density goes as dup/dfl oc out- 
side the inner cutoff, as before. Inside the cutoff, one of 
the distributions (call it Model B-1) has dup/dfl = const 
and the other (call it Model B-2) has dup/dil = 0. These 
models give 



o.7n 



o.ar^SJ-^o 



TTT" 



0.7^ 



U3 
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and 



TTT- 



(12) 



pulsars enclosed within rioff- 

In Model C, we consider a compact population of pulsars 
contained in a solid angle much smaller than any resolu- 
tion solid angle so that dup/dil is effectively a delta func- 
tion. A compact distribution close to Sgr A* cou ld con- 
ceivably arise as a product of dynami cal friction ([Morris! 
[1991 IMiralda-Escude fc Gould [200l . Here we simply 
have 



Nciiy) = Ni 



(13) 



In addition to the above three, one may consider other 
models for the pulsar distribution. For example, the pul- 
sars could be arranged in a central core with a diffuse 
halo. However, most of these other distributions can be 
made as combinations of those we consider. As a result, 
we do not expect the final answers to change by more 
than an order of magnitude. 



within an angular distance oi9 = 25", which corresponds 
to a projected radial distance of r w 1 pc from Sgr A*. 
The allowed ranges for each parameter were chosen to be 
consistent with current measurements and are presented 
below in Table [2| For reference, the best fit values for 
the model with no pulsars present are abh = 0.15 and 
Vf = 0.26 GHz. 

For each grid point in the three-dimensional parameter 
space, we calculate between the measured flux and 
our model and evaluate the likelihood function assuming 
independent Gaussian statistics for measurement errors. 



C{Ni,ahh,Vf 



) = n (2--f ) " ' exp I - jSsM-JoM]' I 



Of cxp [ 



(14) 



The likelihood function is then marginalized over abh 
and Vf to get a distribution for Ni. The marginalized 
likelihood functions for the number of pulsars within 1 
pc of Sgr A* are plotted in Figure [3| 



Table 2 

Searched Parameters for Each Pulsar Distribution Model 




-20 -10 10 20 

Projected Angular Separation from Sgr A* (arcsec) 

Figure 2. Pulsar number density as a function of angular sepa- 
ration from Sgr A* for each of the three distribution models. The 
distributions have been normalized in this figure so that each model 
produces the same number of pu lsars enclosed within the inner par- 
sec (9r ~ 25"). See Section 15^41 for a description of each model. 



5.5. Model Fitting 

To test our distributions, we calculated the values of 
the model given by Equation [H with each of the three 
distributions ag ains t the VLA and G MRT data from 
lAn et all ([200I and lFalcke eFall ([T991 . We fix the pul- 
sar spectral index at apsr = —1.7 and the flux density 
per pulsar at Spsi.uo = 99 fiJy at 1.4 GHz, as described 
in Section [221 The flux from the point source associated 
with the MBH is normalized such that the total measured 
flux of Sgr A* exactly matches the data at 8.45 GHz. The 
spectral index of the MBH point source (abh), the num- 
ber of pulsars (iVi) and the free- free cutoff frequency (;//) 
are allowed to vary. The number of pulsars, A^i, is taken 



Model 


Parameter Ranges" 




Ni 






Uf (GHz) 


A 


[ 0, 50000, 100 ] 


0.05, 0.40, 0.01 




0.05, 1.00. 0.01 




B 


[ 0, 50000, 100] 


0.05, 0.40, 0.01 




0.05, 1.00, 0.01 




C 


[ 0, 5000, 10 ] 


0.05. 0.40. 0.01 




0.05. 1.00. 0.01 





Data given as [ min, max, step size ] 




2000 3000 4000 

Number of Pulsars 



Figure 3. Li kelih ood functions for three pulsar distribution mod- 
els (see Section [5.4l for a description of each model) . The likelihoods 
have been marginalized over the parameters ctbh ^^'^ '^f ^^'^ '1°'^" 
malized so that J C{N) = 1. In each case, the pulsar spectral index 



was taken to be ar. 



-1.7. 



5.6. Results 

The best-fit parameters for the maximum likelihood 
(ML) number of pulsars within the inner parsec of the 
GC for each model is provided in Table [3| and the result- 
ing spectra are plotted in Figure [Jj The number quoted 
in Table [31 is the ML number of pulsars as determined 
from the likelihood distribution and the uncertainties de- 
note the most compact 68% confidence interval around 
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the ML value. If the ML value for a distribution is zero, 
then the upper limits are given at the 68% confidence 
level. In addition to the fiducial pulsar spectral index 
of ttpsr = —1-7, we include spectral indices of —1.0 and 
—2.5. For comparison, the best fit parameters for the 
model with the number of pulsars fixed at zero is also in- 
cluded in Table [3] and odds ratios are calculated against 
this "null" model. 

It is interesting to note that although Models A and B 
give only upper limits to the pulsar population (of ;$10'^ — 
10"), Model C provides a non-zero maximum likelihood 
value of ^10'^. Additionally, Model C provides a better 
fit than the model with no pulsars at all (the "null" model 
in Table O. 

In all of the models considered, an increase in the 
maximum number of pulsars is accompanied by an in- 
crease in the free-free turnover frequency. The increased 
free-free absorption is required to mask the bright low- 
frequency tail of a large pulsar population. Current esti- 
mates of free-free absorption in the region near Sgr A* 
give t urnover frequencies around 330 MHz (jPedlar et alJ 
Il989t) . Estimates of the free-free absorption of the flux 
from Sgr A*, however, assume a power-law fiux density 
for Sgr A* and would likely underestimate the absorp- 
tion if a pulsar population were present. As a result, 
any independent measurement of the free- free absorption 
along the line of sight of Sgr A* that does not assume a 
spectrum of the Sgr A* source could place an important 
constraint on the pulsar population in the inner parsecs 
of the Galaxy. 

From current radio measurements of the inner parsec of 
the GC, total pulsar populations (that is, both CPs and 
MSPs) of up to ^10^ are consistent with observations, 
regardless of the underlying spatial distribution. 

Table 3 

Number of Pulsars within 1 pc for Given Model and Spectral Index 



Model Parameters" 



Model 


Qpsr 






(GHz) 


Ndof 




Odds 


A 


-1.0 
-1.7 
-2.5 


< 16.2 

< 7.4 

< 2.2 


0.14 
0.14 
0.14 


0.37 
0.38 
0.38 


10 
10 
10 


1.08 
1.16 
1.17 


10-0.64 

10-1.05 
10-1.35 


B-0 


-1.0 
-1.7 
-2.5 


r, Q+2.5 

< 4.9 

< 0.9 


0.17 
0.15 
0.13 


0.45 
0.53 
0.47 


10 
10 
10 


0.79 
1.07 
1.27 


10-0.55 
10-1.21 
10-1.42 


B-1 


-1.0 
-1.7 
-2.5 


< 3.8 

< 1.5 

< 0.4 


0.14 
0.13 
0.14 


0.40 
0.41 
0.38 


10 
10 
10 


1.06 
1.14 
1.17 


10-0.85 
10-1.03 
10-1.22 


B-2 


-1.0 
-1.7 
-2.5 


< 2.4 

< 1.2 

< 0.4 


0.14 
0.14 
0.14 


0.36 
0.38 
0.38 


10 
10 
10 


1.07 
1.13 
1.17 


10-1.02 
10-1.00 
10-1.23 


C 


-1.0 
-1.7 
-2.5 


1.1 ±0.4 

-'-■-'--0.6 


0.25 
0.21 
0.17 


0.40 
0.47 
0.58 


10 
10 
10 


0.55 
0.59 
0.75 


10+0.28 
10+0.12 
10-0.46 


Null' 







0.15 


0.26 


11 


0.81 


1.0 



Best fit model parameters for maximum likelihood (ML) number of pulsars. 
If the ML number of pulsars is zero, then the 68% confidence upper limit is 
used and reported with a "<". 
^ The "Null" case fixes the number of pulsars at zero. 



Projected Radius of Effective Resolution at Sgr A* (pc) 
^ ^ 1.0 0.3 0.1 0,03 0.01 




0.1 0.3 1.0 3.0 10.0 



Frequency (GHz) 



Figure 4. Top: Fit curves for model pulsar dis tributions to the 
VLA (squares) and GMRT (diamo nd) data from I An et al.l (|2005l) 
and VLA (filled circles) data from lFalcke et all 1)19981 ). Since the 
best fit for Models A and B-1 indicate that the most probable 
number of pulsars is zero, we have instead plotted the best fit 
assuming the 68% confidence upper limit number of pulsars are 
present. The upper axis gives the projected radial distance of the 
effective resolution bea m a t the distance of Sgr A* using the VLA 
(for details, see Section fS.SH . Note that although the fits were made 
using both the VLA and GMRT data, the curves above only apply 
to the VLA data points (squares and filled circles). Middle and 
Bottom: Components to the observed fiux of Sgr A* from the MBH 
point source and surrounding pulsars for Models A and C. Note 
how the pulsar component in Model A rises quickly with decreasing 
frequency as a result of both the pulsar spectrum and the increasing 
number of pulsars in the beam. The total unabsorbod flux is also 
shown for comparison. 



6. DIFFUSE GAMMA-RAY EMISSION AND MSPS 

MSPs arc known gamma-ray sources ()Abdo et al.l 
l2009a[) . As a result, we may set constraints on the MSP 
population in the GC by measuring the diffuse gamma- 
ray emission from the GC. In a recent analysis of the 
first two ye ars of data from the Fermi Gam ma- Ray Space 
Telescope, iHooper fc Goodenoughl ()201lL herafter, HG) 
observed an excess of gamma-ray flux toward the inner 1° 
(150 pc) of the GC, with a significant excess within 0.25° 
(w 40 pc). HG argued that the signal is consistent with 
the annihilation of 7-10 GcV dark matter particles with 
a cusped halo distribution around Sgr A* a nd difficult 
to exp lain using known astrophysical sources. lAbazaiianl 
(|2011f ). however, claimed that while the spectrum may be 
inconsistent with the average pulsar spectrum, it is con- 
sistent with some pulsar spectra and therefore could be 
explained by astrophysical sources. HG provide a spec- 
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trum for the gamma-ray excess in the GC of the form 



cx E-^ exp{-E/Ec 



(15) 



0.99;°;J° and E.^t 



1.92t°;2i GeV. These 



where T 

parameter values are consistent (to current uncertain- 
ties) with the spectra of 16 out of 46 pulsar s in the first 
Fermi LAT catalog of gamma-ray pulsars (|Abdo et al.l 
l2010bD . We proceed assuming that the observed excess 
seen by HG is real (however, see iBovarskv eralllMTI ). 
follows the spectrum fit by HG and is entirely caused by 
a collection of MSPs in the GC. 

The spectrum is normalized so that E'^dN^/dE « 
10"^ GeV cm-2 s'^ at E = 1.0 GeV based on the HG 
plots, giving 



dN^ 
'dE 



= 1.7 X 10" 



X exp 



GeV"^cm"2g-i 



E 



E 



1 GeV 



-0.99 



1.92 GeV 



(16) 



Calculating the integrated energy flux from the above 
spectrum over 0.1-100 GeV as 



S-y • — - 



100 GcV 



E- 



0.1 GcV 



^dN^ 
'dE 



-dE 



(17) 



and letting = Aird'^ fnS^, we find that the gamma- 
ray luminosity of a source at the GC is « 4 x 
Iq36 Qj-g s~^/q. The correction factor, fa, is similar to 
the beaming fraction in radio pulsars a nd is generally 
taken to be unity in most modern models (jWatters et al.l 
[20091) . 

Following calculations mad e to estimate the n umber of 
MSPs in globular clusters (|Abdo et al.l l2010a|) . we can 
estimate the MSP population of the GC as 



.^MSP 



(18) 



where E is the total spindown luminosity of a MSP, 
rjy = Ly/E is the "efficiency" of converting spin-down 
power into gamma-rays, and angled brackets denote an 
average over the population. The values of (E) and {r]-y) 
are taken from the local {d < 1 kpc) MSP population 
to avoid selection effects. The mean spin-down luminos- 
ity is taken to be (E) = 1.1 x lO'^"^ erg s~^ from the 27 
MSP s within 1 kpc of Earth listed in the ATNF cat- 
alog (jManchester et al.]l2005f ). Since (77.^) oc oc d^, 
the uncertainties in {ij^) are dominated by distance un- 
certainties and values range from (77.^) = 10~'^ — 1 in 
the fir st Fermi gamma-ray pulsar catalog (|Abdo et al.l 
l2010b[) . To mitigate this distance problem, we take only 
those seven nearby pulsars in the catalog for which dis- 
tances could be measured accurately with parallax. For 
these pulsars, we find (77^) = 0.08 ± 0.04. Using Equa- 
tion [181 we get an estimate of A'msp ^ 5000. In order to 
contribute the observed excess, these pulsars would be 
located within 1° (150 pc) of Sgr A*, with the highest 
concentration within 0.25° (40 pc). 

Our estimate for the number of MSPs in the GC is 
essentially an upper bound for the population. How- 



ever, modeling the background component of the diffuse 
gam ma-ray emission in the GC is still somewhat uncer- 
tain (|Abdo et al.ll2009bD and, as a result of this uncer- 
tainty, our estimate can only be taken as an approximate 
upper limit. Additionally, as this is a measure of the 
excess gamma-ray flux in the region (with the Galactic 
plane and a central point source coincident with Sgr A* 
subtracted), there exists the possiblity that a signicant 
number of MSPs are unaccounted for in this estimate. 
Assuming the HG gamma-ray excess in the GC is real 
and does not suffer from systematic errors in background 
subtraction, we see that it is not inconsistent with a cen- 
trally concentrated population of -^10^ MSPs in the inner 
tens of parsecs from Sgr A*. 

Finally, we note that similar estimates for th e MSP popu- 
lation in the GC have been made in the past. I Wang et al.l 
()2005D used a model proposed bv lZhang fc Chengj (|1997[ f 
to predict the emission of gamma-rays from MSPs as a 
function of global pulsar parameters like spin period and 
polar magnetic fleld. They used Monte Carlo methods 
to simulate a population of pulsars and measured the to- 
tal gamma-ray luminosity. Com paring this luminos ity to 
measured values from EGRET. IWang et all ()2005D esti- 



mated a GC MSP population of A'^msp 
EGRET field (r « 1.5°, 220 pc). 



6000 in the 



7. MASSIVE STARS IN THE GALACTIC CENTER 

The central parsec of the GC is one of the most active 
massive star formation regions in the Milky Way and 
is cur rently known to co ntain about 200 young massive 
stars (IGenzel et al.l[2?)10h . As massive stars are the pro- 
genitors of neutron stars (NSs), we may use current stel- 
lar populations to estimate the number of pulsars in this 
region. 

The inner parsec stellar population is divided into three 
fairly distinct regions. The innermost region {R < 1") is 
the so-called "S-star Cluster" of main sequence B-stars. 
These stars can be fitted with a standard Salpeter initial 
mass function (IMF) for a single star formation event or 
for a continuous st ar-forming populati on with ages of a 
few Myr to 60 Myr (jBartko et al.ll2010D . Outside this re- 
gion (1" < R < 12"), the stars are largely arranged in at 
least one disk (possibly two) of mass M ^ 10"^ Mq with a 
top-heavy IMF of dN/dm oc m,-^-'^^ (jBartko et al.ll2010( ). 
The early-type stars in this region appear to have been 
formed in a starburst ^6 Myr ago. Outside the disk re- 
gion (R > 12") the stellar population is again consistent 
with a Salpeter IMF. 

Foll owing similar calculations bvlLazio fc CordesI ()2008[ ) 
and iFaucher-Giguere fc Loebl (|2011[ ). the current stellar 
populations can be used to estimate the number of ra- 
dio pulsars harbored in the central parsec. The number 
of active CPs and MSPs beamed toward Earth are esti- 
mated to be 



CP 



fpsrfbfrfv-^n 



and 



A^MSP " fpsrfbfrfvfrNn 



(19) 



(20) 



respectively, where iV„s is the number of NSs in the inner 
parsec, fp^r is the fraction of NSs that form pulsars, /t, = 
0.2 is the beaming fraction, /t is the fraction of pulsars 



Constraints on GC Pulsar Populations 



11 



with ages less than the typical pulsar radio lifetime, /„ is 
the fraction of pulsars with birth velocities small enough 
to be retained in the inner parsec and fr is the fraction 
of NSs that are recycled into MSPs. The fraction of 
NSs that form pulsars is taken t o be f^.^r ^ li althoug h 
this factor is still fairly uncertain (jLazio fc Corde£ll2008 ). 
The terms Nns, /t, fv and fr are discussed below. 

7.1. Neutron Star Population in the Inner Parsec (N„s) 

The total number of NSs residing in the inner parsec 
may be estimated from current observations of the mas- 
sive star population. Velocity measurements of stars al- 
low for the total dynamical mass enclosed within a few 
parsecs of Sgr A* to be determined. Subtracting the 
mass of the MBH gives an extended mass of ^10^ M, 
in stars within a p arsec of Sgr A* (|Genzel et al.l[l9& 
iSchodel et aT1[2009[ ). With the total mass of stars estab 
lished, the number of NSs can be calculated for a given 
initial mass function (IMF) and mass range of stars that 
end their lives as NSs. 

A typi cal NS progenitor mass range is 9Mq < M < 
25Mq (|Heger et al. 2005). However, such ranges are the- 
oretically determined for only isolated non-rotating stars 
and rely on mass-loss and stellar wi nd models that ar e 
poorly constrained by observations (jHeger et al.l 120031 ). 
If very massive stars have higher mass-loss rates than 
modeled, then NSs could form from stars with initial 
masses above > 25Mq. That this may be the case is 
supporte d by the lim i ted ob servational constraints. For 
example. iMuno et al.l (|2006D detected an X-ray pulsar in 
the young Galactic cluster Westerlund 1 that requires a 
progenitor mass of M > AOMq to have formed in the age 
of the cluster. 

To calculate the number of NSs produced by the GC 
massive star population, we consider IMFs of the form 
dN/dm oc m~" over a range of masses from O.IMq 
to IOOMq. For the case of the standard Salpeter IMF 
(a = 2.35) observed in most of the inner parsec and 
a progenitor mass range of 9Mq < M < 25Mq, a to- 
tal of Nns ~ 4900 NSs are produced. If the top-heavy 
IMF {a = 0.45) observed in the young disks is taken to 
hold for the whole central parsec, a total of Nns ~ 5700 
are produced. If a wider range of progenit or masses is 
taken, say 9Mq < M < 40Mq to allow the iMuno et all 
()2006D observation, the total NS populations increase to 
Nns ~ 5700 and Nns ~ 9500 for the Salpeter and top- 
heavy IMFs, respectively. Even given a wide range in 
IMF and progenitor mass, the above results are consis- 
tent to an order of magnitude with a NS population of 
Nns ~ 10^. In terms of NSs formed per stellar mass, we 



find a value of /3„ 



lO-^M-i for the -lO^Mo worth 



of stars within a parsec of Sgr A*. 

7.2. Fraction of Still-Active Pulsars (fr) 

The typical lifetimes over which pulsars maintain ac- 
tive radio emission are r ~ lO'' yr for CPs and r ~ 
109 - iQio yr for MSPs. As a result, the fraction of 
pulsars formed recently enough to still be active must 
be considered in estimates of the observable population. 
This factor will depend on the star formation history 
(SFH) of the region. In the case of continuous star for- 
mation, the active fraction of pulsars can be estimated to 



be fr ~ r/tsf, where tgf is the amount of time elapsed 
since the star formation began. 

7.3. Fraction of Pulsars Retained (fy) 

From the high observed velocities (~10^ — 10'^ km s^^) of 
some pulsars, it has been inferred that NSs are given a 
large "kick" velocity at birth as a result of binary dis- 
association or an asymmetric supernova explosion (or 
both). A pulsar created in the inner parsec will be re- 
tained in the inner parsec only if its birth velocity does 
not exceed the local escape velocity of its orbit around 
Sgr A*. Assuming only the influence of a MBH of mass 
4 X 10^ Mq, the escape velocity at a distance r 
from Sgr A* is given by 



185 km s 



M 



4 X 106 M, 



1/2 



1 pc 



-1/2 



(21) 

Assuming a Maxwellian distribution of birth vc- 
locities with a mean of (vh irth) — 380 km s~^ 
(jFaucher-Giguere fc Kaspill2006l) . we find retention frac- 
tions of fv ~ 0.1 and fv ~ 0.25 for radial distances of 
r = 1 pc and r = 0.5 pc, respectively. However, the ac- 
tual shape of the pulsar birth velocity distribution is not 
well constrained and one may worry that the Maxwellian 
dist ribution is arbitrary. In an analy sis of pulsar veloci- 
ties, iFaucher-Giguere fc Kaspl (|2006f ) consider six differ- 
ent pulsar velocity distributions. Repeating our calcu- 
lation for each distribution, we find retention fractions 
in the ranges of fv ~ 0.05 — 0.4 and fv « 0.1 — 0.5 for 
distances of r = 1 pc and r = 0.5 pc, respectively. 
If no other effects are important, we would expect a re- 
tention fraction of /.„ 0.1. However, a similar analy- 
sis to the one performed above would underestimate the 
retained NS populations in globular clusters by orders 
of magnitude. Observations of pulsars and low-mass X- 
ray binaries (LMXBs) suggest that up to ^^10% of all 
NSs formed in so me globular clusters may be retained 
(jPfahl et al.l(200l . However, current models for isolated 
pulsars predict that < 1% of NSs will have birth ve- 
locities below the < 50 km s~^ globular cluster escape 
velocities. This "retention problem" is currently an un- 
solved problem, but likely has to do with the high stellar 
densities and bi nary fractions fou nd in the cores of globu- 
lar clusters fsee lPfahl et al.ll2002l and references within). 
Since the GC has even higher stellar densities than cores 
of globular clusters, we expect a similar heightening of 
the retention fraction and thus adopt as a nominal value 

7.4. Fraction of Pulsars Recycled to MSPs (fr) 

MSPs are thought to be formed when a NS in a binary 
gains angular momentum through a ccretion of matter 
in a process known as "recycling" (lAlpar et al.l 119821 : 



iBhattacharva fc van den Heuvellll99H) . Thus, the frac 
tion of NSs recycled into MSPs must be determined to es- 
timate the MSP population. Within 3 kpc of the Sun, the 
birthrate of pulsars with 400 MHz pseudo-luminosities 
above 1 mJy kpc^ is observed to be ^10~ ^ yr~^ and 
>10~6 yr-i for CPs and MSPs, respectively (|Lvne et al.l 
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Il998f ). Using these birthrates, wc may infer that the re- 
cychng fraction is at least /,. > 10"'^ in the Galactic disk. 
The increased stellar density and stellar encounter rate 
in the GC will very likely increase this fraction. LMXBs, 
the assumed progenitors of MSPs, have been found to 
be ~100 times more abun dant in globular clusters than 
the general Galactic field (jGlarkl 119751 : iKat^ 119751 ). As 
the central parsec of the GC has a higher stellar density 
than globular clusters, we would expect at least a simi- 
lar overabundance of LMXBs and their resultant MSPs 
as is seen in globular clusters. As a result, we adopt a 
recycling fraction of fr ^ 0.1. 

7.5. Pulsar Estimates for Various Star Formation 
Histories 

As the star formation history (SFH) of the central par- 
sec of the GC is still somewhat uncertain, we will con- 
sider two general SFHs suggested by current observa- 
tions. In the first case, we take the massive-star disk(s) 
to have formed in a well-defined starburst ~6 Myr ago 
(jPaumard et all 120061: iBartko et all I2O1O0 and assume 
that the rest of the central parsec has experienced con- 
tinuous star formation over the age of the Galaxy. In 
the second case, we consider SFHs based on spectro- 
photometry of cool giant stars which indicate that most 
of the stars in the inner parsec were formed > 5 Gyr ago 
but also sh ow an increased star formation rate i n the last 
~100 Myr (|Blum et al.ll200a IPfuhl et al.llMl . 

7.5.1. Continuous Star Formation + Disk Starburst 

In the first SFH, pulsars can come from both the general 
population of stars in the central parsec and the young 
disk population. For the general population, we take the 
total mass of stars to be M ~ lO^M© and assume con- 
tinuous star formation over the last ^10^° yr. From the 
parameters discussed above, the total number of NSs is 
found to be iV„, - {^n^ / W^M'^) . Taking a CP 
active radio lifetime of r 10^ yr, continuous star for- 
mation over ^10^" yr will give the fraction of CPs still 
active to be fr ^ 10^^. From Equation [TOl the CP con- 
tribution from the continuous star forming region of the 
inner parsec is 



NET 




(22) 



For MSPs with a radio lifetime of ^10^° yr, the fraction 
still active is 1. Adopting the recycling fraction 

of fr ^ 0.1 discussed above, the contribution of MSPs 
from the general population of the inner parsec is given 
by Equation [20l as 



/vs<=" ^200 



Additionally, the contribution of pulsars from the 
massive-star disk(s) must be considered. The disk is as- 
sumed to have a population of stars formed in a starburst 
event ^^6 x 10^ yr ago with a total stellar m ass of ~1O^M0 
(jPaumard et al.ll2006l: feartko et al.ll2010D . Since the age 
of the disk is comparable to the active radio lifetime of a 
CP, the fraction of CPs still active is taken to be /r ~ 1. 




From Equation [THl the CP contribution from the disk is 
found to be 



Ardisk 




isk 



\ 104 Mr. 



(24) 



The disk population is not expected to produce any 
currently observable MSPs as the short timescale of 
^6 X 10^ yr provides insufficient time to create and evolve 
a NS population into MSPs. As a result, this first SFH 
produces roughly Nqp ^ 20 CPs and TVmsp ^ 200 MSPs. 

7.5.2. SFH from Observations of Cool Ciant Stars 

The SFH has also been estimated by comparing simu- 
lated populations with the observed cool giant stars in 
the central parsec. Such simulations allow the average 
star formation rate to be calculated as a function of 
look-back time for a few coarse t ime b ins. In two sep- 
arate analyses, both lBlum et al.l (|2003[ ) and IPfuhl et al.l 
(|20Tl[) found that >80% of the stehar mass in the cen- 
tral parsec was formed >5 Gyr ago and that there has 
been increased star f ormation in t he las t ~100 Myr. In 
their best fit models, I Blum "eTall (pool found an aver- 
age star formation rate of ~3 x 10~^ M(^yv~^ with i n 2 pc 
of Sgr A* from 10 to 100 Myr ago and IPfuhl et"aI1 ()2011l ) 
found an average star formation rate of ~lO~'^M0yr~^ 
within 1 pc of Sgr A* from 50 to 200 Myr ago. Both cases 
are consistent with ^IO^AIq worth of stars being formed 
in the inner parsec of the GC in the last '^100 Myr. 
If there has been continuous star formation in the last 
~100 Myr, then the fraction of still active CPs would 
be fr ^ 0.1. Taking all other parameters as before, the 
continuous formation of ^IO^A/q worth of stars over the 
last ~10^ yr would produce 



^cp - 20 ( 



' M{t< 10* yr)~ 



0.01 M 







105 Mp 







(25) 



If the recent star formation all took place in the last 
^10^ yr, then the fraction of CPs still active would be 
fr ^ 1- In this case, the number of active CPs would be 




Arburst 
^^CP 



In either of the above cases, the majority (>90%) of the 
star formation took place at look-back times >10* yr 
ago. As a result, the number of MSPs produced will 
be approximately the same as the first SFH considered, 
namely iVMSP ^ 200. 

7.5.3. Upper Limits to the Pulsar Population 

Using the above estimates of CP and MSP populations 
for a range of observationally supported SFHs, upper lim- 
its may be set on the total allowable number of pulsars in 
the inner parsec. For CPs, the most favorable formation 
scenarios produce Nqp ~ 200. For MSPs, a variety of 
SFHs consistently produce TVmsp ~ 200 (/,-/0.1). Since 
the recycling fraction is unknown for the extreme condi- 
tions of the inner GC, an upper limit of TVmsp ~ 2000 
may be set by adopting fr ^ 1. Thus, observations 
of current stellar populations place an upper limit of 
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a few X 10^ on the number of active radio pulsars beamed 
toward Earth in the inner parsee of the GC. 

8. PULSAR WIND NEBULAE IN INNER 20 PC 

Using a total of 1 Ms of Ch andra ACIS- I obser vations of 
the inner parsecs of the GC. lMuno et"all (pOOSi ) eompiled 
a catalog of 34 diffuse X-ray emitting features. Based on 
the X-r ay luminosities and sizes of the sources in their 
catalog. [Muno et al.l (I2008D expect ~20 pulsar wind neb- 
ulae (PWNe) to be present within 20 pc of Sgr A*. Since 
PWNe are powered by pulsars, we may use this inferred 
population of PWNe to estimate the pulsar population 
in the inner 20 pc of the GC. 

Pulsars can lose their rotational kinetic energy by the 
release of relativistic winds of charged particles. The 
winds exert a pressure upon and deposit energy into 
the surrounding interstellar medium, producing lumi- 
nous PWNe th at radiate acro ss the electromagnetic spec- 
trum (see. e.g.. IGaensler fc Slang. 20061 ) . As a result, the 
luminosity of the PWN will be directly related to the 
spin-down luminosity of the pulsar given by 



E = -dErot/dt = 47r2/P/p3, 



(27) 



where / and P are the moment of inertia and period, 
respectively, of the pulsar. 

Since the spin-down luminosity of a pulsar decreases 
with increasing age, one would expect the most luminous 
PWNe to contain young pulsars (however, older pulsars 
may be "r ecycled" sufficient ly to power PWNe as de- 
scribed bv lCheng et al.|[2"006( ). Of the 30 confirmed pul- 
sars associated with PWNe, 25 have characteristic ages 
(tc = P/2P) of Tc < 10^ yr and 23 have characteris- 



tic ages Tc < 10'^ yr (iRobertsI [20041 ) . Thus, we adopt a 



pwn 



10" 



yr- 



lifetime for a typical PWN of t 
Given the observed number of PWNe in the GC {Nobs ^ 
20) and a typical lifetime of tp^n ^ 10^ yr, we find a 
mean rate of formation of PWNe over the last ^10^ yr 
to be 



pwn 2 X 10 yr 



20 



10^ yr 



(28) 



Assuming that the PWN formation rate has remained 
constant over the last ~10^ yr, we may estimate the num- 
ber of CPs in the inner 20 pc to be 



Nqp ~ PpwnTpsr fbfv fpwm 



(29) 



where r^sr ^ 10^ yr is the typical radio lifetime of a CP, 
fb = 0.2 is the beaming fraction, is the fraction of 
pulsars with birth velocities low enough to be retained 
in the inner 20 pc, and fpwn 1 is the fraction of pulsars 
that form PWNe. 

Taking t he mass of the centr al 20 pc to be M = 3 x 
10^ Mq (jLindgvist et al.|[l992[ ). we find that pulsars must 
have velocities vurth < 115 km to remain gravita- 
tionally bound to the inner 20 pc. From the birth veloc- 
ity dis tributions considered bv lFaucher-Giguere fc Kaspil 
(|200a . we find /„ w 0.05 - 0.30. However, the distance 
traveled by a pulsar is given by 



10 pc 



100 km s" 



105 yr 



(30) 



Thus, pulsars with velocities high enough to become 
gravitationally unbound will also have velocities high 
enough to escape the inner 20 pc on timescales com- 
parable to the PWN Ufetime. As a result, the PWNe 
observed within the inner 20 pc are very likely to remain 
there and we can take fy ^ I. 

From Equation [29l we see that if the PWN birth rate has 
been constant over the last ~^10^ yr. we would expect 
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(31) 



CPs within 20 pc of Sgr A*. 

Finally, we note that of the ^20 PWN candidates iden- 
tified bylM uno et all (|2008[ ). 4 fall within a projected ra- 
dial distance of 1 pc from Sgr A*. Assuming the number 
of pulsars scales accordingly, then Equation [31] predicts 
Ncp ~ 80 CPs within the inner parsec of the GC. 

9. SUPERNOVA RATE IN THE GALACTIC CENTER 

Neutron stars are formed as the end products of core- 
collapse supernovae (CCSN). An estimate of the rate of 
CCSN in the GC would therefore offer a constraint on the 
pulsar population. The CCSN rate is estimated below for 
both r < 150 pc and r < 20 pc. 

9.1. CCSN Rate Within r < 150 pc of Sgr A * 

By measur i ng th e total mass of ^^Al in the Galaxy, 
iDiehl et all ()2006D estimate the Galactic CCSN rate to 
be PccSN = 1-9 ± 1.1 century"^. One may, in princi- 
ple, scale this estimate to smaller regions of the Galaxy 
using massive star populations. Taking the inner 500 
pc t o contain 10 % of the Galaxy's massive star forma- 
tion (jFigedl20"08[ ) we can estimate that the inner ^--^150 pc 
contains ~2% of the massive star formation and therefore 
should have a CCSN rate of Pccsn ~ 0.04 ccntury^^. 
iCrocker et al.l (|2011[ ) estimate a similar rate and show 
that it is consistent with SN rate estimates from in- 
frared observations, stellar composition. X-ray emission, 
gas turbulence, an d high-velocity compact clouds (see 
ICrocker et al.l[20Tll . and references within) . We may now 
estimate the CP population in the GC to be 



A'^CP — fpsrfbfvTpsrPcCSN 



(32) 



where Pccsn « 4 x 10 ''yr ^ is the CCSN rate, Tpsr ^ 
10^ yr is the mean canonical pulsar lifetime, fb — 0.2 
is the fraction of pulsars beamed toward Earth, fy is 
the fraction of pulsars with birth velocities small enough 
to be retained by the GC and fpsr 1 is the frac- 
tion of CCSN tha t result in active pulsars. Using the 
distributions from iFaucher-Giguere fc KasDil (I200I . the 
fraction of pulsars with birth velocities smaller than 
the escape velocity Ve ~ 200 km s^^ at 150 pc ranges 
from fy w 0.1 — 0.4. These values give an estimate of 
iVcp '^100. Likewise, accounting for the longer ages for 
MSPs {Tpsr ^ 10^° yr), we can estimate the MSP pop- 
ulation to be A^MSP ~ 10^ /r, where fr is the fraction of 
NSs that get recycled to MSPs (see Section [73) • 

9.2. CCSN Rate in Inner 20 pc from X-ray 
Observations 
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Studies of diffuse X-ray emission can also provide in- 
sight into the SN rate in the GC. Using over 60 ks of 
Chandra ACIS-I observations. iMuno et all (I2004D found 
that the diffuse X-ray emissions in the GC could be 
explained by a two-temperature plasma composed of a 
"soft" component {kT k, 0.8 kcV) and a "hard" compo- 
nent {kT K, 8 keV). Assuming the soft component of 
the plasma is primarily heated by SNe, an estimate for 
the SN rate can be made by observing the loss of energy 
from the inner 20 pc. 

Let us first consider the case in which the soft component 
of the plasma just cools radiatively. The X-ray luminos- 
ity of the soft component of the plasma in th e inner 20 pc 
is Lx « 3 X 10^6 erg s'^ (|Muno et al.ll2004n . If each SN 
transfers ^1% of its total kinetic energy of ^10^^ erg 
to the plasma, then a SN rate of (3sn ~ 10^^ yr~^ 
is required to maintain the currently observed temper- 
ature. Taking this rate to be constant and /„ ~ 0.1 (see 
Section ^ gives an estimate for the CP population of 
A^cp ^2(/„/0.1). 

If the plasma is unconfined, it can also cool through adia- 
batic expansion. Rough estimates put this cooling rate at 
Larf « 9 X 10^8 erg s'^ (jMuno et al.ll2004r ). which would 
require a SN rate of Psn ~ 3 x 10~^ yr~^. Again as- 
suming this rate is constant and /„ ^ O.l as above, the 
estimated CP population is A^cp ^ 600 (/t,/0.1). 
The above estimates assume that the SN rate is con- 
stant over the radio lifetime of a CP (-^10^ yr). How- 
ever, the SN rate estimates are only required to hold 
over a characteristic cooling time of ic E/L^ where 
E ^ b X 10^° erg s~^ is the total thermal energy stored 
in t he plasma and L i s the appropriate cooling luminos- 
ity ()Muno et al.ll200^ . The cooling timescales for radia- 
tive cooling and adiabatic expansion are 5 x 10^ yr and 
2 X lO'' yr, respectively. 

Keeping this caveat in mind, the assumption of a con- 
stant SN rate does allow useful bounds to be put on the 
pulsar populations in the GC. Taking a constant SN rate 
from radiative cooling to be a lower bound and a con- 
stant SN rate from adiabatic expansion to be an upper 
bound, we find the CP population in the inner 20 pc to be 
1 ^ ^CP < 10^. Likewise, for MSPs with t - 10^° yr, 
we estimate a population of IQ^fr S ^msp ^ 10® /r, 
where fr is the fraction of NSs recycled into MSPs (see 
Section [73). 

10. DISCUSSION 

10.1. Summary of Estimates 

We have used observations over a wide range of wave- 
lengths to make order of magnitude estimates of the num- 
ber of pulsars allowed within r < 1 pc and r < 150 pc of 
Sgr A*. The estimates are summarized in Tabled 

10.1.1. Pulsar Population Withm r < 150 pc 

Limits on the pulsar population within 150 pc of Sgr A* 
may be set using the five known pulsars in the inner 15', 
the ca talog of compact radio sources by iLazio fc CordesI 
(l2008i) . measurements of an ex cess gamma-ray flux 
bv iHooper fc Goodenoughl (|20Tll) and estimates of the 
Galactic core-collapse supernova rate. 
The five known pulsars in the inner 15' provide the 



strongest current evidence for an intrinsic pulsar pop- 
ulation in the GC region. A M onte Carlo population 
analysis by iDeneva et al.l (I2009D showed t hat t he pulsar 
detections in the survey bv IDeneva et al.l (|2009[ ) indicate 
a population of at least N > 100 pulsars within 100 pc 
of Sgr A*. 

The catalog of comp act radio sources compiled by 
ILazio fc CordesI (pOOSi ) allows for upper limits to be 
placed on the pulsar population in the inner 150 pc. 
The VLA survey produced a total of 170 compact steep- 
spectrum sources. Although pulsars can not be unequiv- 
ocally classified in an imaging survey, ILazio fc CordesI 
(|2008| ) estimate that A^ 10 of the sources were likely 
pulsars. Using this estimate and the survey sensitivities, 
a conservative upper bound of N < 10^ may be set on 
the pulsar population in this region. Even if it turns 
out that there are no pulsars in the catalog, an upper 
limit oi N < 10^ may be set to 99% confidence level. We 
note that although these upper limits should hold for the 
entire survey region, the increasing background temper- 
ature and sidelobes caused by extended emission toward 
the inner GC mean that the inner field of the survey 
(half-power radius of 15') likely experienced decreased 
sensitivity and could potentially hide a significant pulsar 
population in the immediate vicinity of Sgr A*. 
Using the first two ye ars of Fermi data, 
IHooper fc Goodenoughl (|2011| ) claim to have detected 
an excess diffuse gamma-ray flux in the inner 150 pc 
of the Galaxy, which they attribute to annihilating 
dark matter particles. If we assume the gamma-rays 
come instead from a collection of MSPs, an upper 
limit to the number of MSPs in the GC may be set. 
Following s i milar calculations for globular clusters by 
lAbdo et al.l(l2010aD . we find that the excess is consistent 
with a population of ^^5 x 10^ MSPs. This estimate is 
nominally an upper limit, but considering the systematic 
uncertainties in current GC gamma-ray background 
models and difficulties subtracting the point source 
associated with Sgr A*, we adopt this value as only a 
lower limit on the upper bound of MSPs in the inner 
150 pc. 

As radio pulsars are formed in corc-coUapse SNe, the 
SN rate also provides a constraint on the pulsar popu- 
lation in the GC. By scaling down the Galactic SN rate 
based on r nassive star populatio ns in a manner similar 
to that of iCrocker et al.l ()201lD , we estimate a SN of 
PccSN ~ 0.04 century"^. Such a SN rate indicates a 
CP population of A^'cp ^ 100 and an MSP population 
of A^MSP ^ where fr is the recycling fraction dis- 

cussed in Section [721 Since the recycling fraction of field 
pulsars is fr > 10"'^ and potentially as high as fr 0.1 in 
globular clusters, we adopt a nominal value oi fr ~ 10~^, 

so iVMSP ~ 103(/,/10-2). 

Finally, we note that all of the upper limits consid- 
ered for pulsar populations on the large scale of 150 pc 
may not include contributions from the few inner par- 
sees as a result of decreased sensitivity or failure to 
incorporate recent starburst activity. As a result, we 
consider these limits mainly applicable in the region 
^ T>c r < 150 pc. The above estimates are consis- 
tent with a total pulsar population (that is, both CPs 
and MSPs) of 10^ < Npsr < 10^, a CP population of 
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Table 4 

Summary of Pulsar Population Estimates by Method 







r < 1 pc 






r < 150 pc 


Method 


CP 


MSP 


Total 


CP 


MSP Total 



Pulsar Surveys 






>102 


Radio Point Sources 






<10* 


Radio Spectrum A 




< 7.4 X 10-' 




Radio Spectrum B-1 




< 1.5 X 10-^ 




Radio Spectrum C 




(1.1 ±0.4) X 10=* 




Diffuse Gamma-ray 






~5 X 10^ 


Massive Stars" 


<102 <10^ 






PWNe 


10^ 






Supernovae"' 


<10=* <10^(/r/10-2) 


<10'' 


10^ 10^ (/r/10~2) 



^ Typical values for the "recycling fraction" are fr > 10 in the Galactic field and potentially as 

high as fr 0.1 in the central parsec (see Section 17.411 
^ Supernovae estimate given at r < 20 pc and r < 150 pc 



^CP < 10^ and an MSP population of Nmsp < 10^. 

10.1.2. Pulsar Population Within r < 1 pc 

Limits on the pulsar population in the inner parsec of the 
GC have been set using the non-detections of high fre- 
quency directed pulsar searches, the spectrum of Sgr A* 
on arcsccond scales, the population and star-formation 
history of the massive star progenitors of NSs, the obser- 
vations of ~20 PWN candidates in the inner 20 pc, and 
the limits on the SN rate based on X-ray observations in 
the inner 20 pc. 

Directed pulsar searches of the inner parsec of the Galaxy 
have been conduct ed with the GBT a t frequencies of 
5, 9 and 15 GHz (jPeneva et all 120091 : iMacquart eraD 
mmh . Since none of these searches made any detections, 
an upper limit may be set on the total pulsar popula- 
tion that is still consistent with a null result. Us i ng the 
surv ey parameters p r ovide d in the iDeneva et al.l (|2009f ) 
and IMacquart et aLl (|201Clf ) surveys, we find that up to 
Npsr ^ 10'^ pulsars (both CPs and MSPs) may be present 
in the inner parsec. 

Since the compact radio source Sgr A* is broadened by 
interstellar scattering («1" at 1 GHz), the observed flux 
density may actually be a combination of the emission 
near the MBH and a diffuse component from a popula- 
tion of unresolved pulsars. By requiring that this two- 
component system reproduce the observed spectrum of 
Sgr A*, constraints may be placed on the pulsar popula- 
tion on arcsecond scales. We consider a variety of spatial 
distributions and find that a total population of '^10'^ 
pulsars is consistent with flux density measurements, re- 
gardless of spatial distribution. The existence of such a 
large pulsar population would distort the low-frequency 
measurements of the intrinsic spectrum of Sgr A* and 
the free- free absorption along the line of sight of Sgr A*. 
Upper limits to the total pulsar population can also be 
set by studying the populations of massive stars that end 
their lives as NSs. Infrared observations of the present 
day population of massive stars allows for estimates of 
the star formation history. We consider the case of two 
general SFHs and find that the CP population can get as 
high as A^cp ^ 200 only under the most favorable con- 
ditions. More typical estimates for the CP population 
are A'cp ~ 20, with most of these being formed in the 



young disk of massive stars located «0.5 pc from Sgr A*. 
The MSPs arc less sensitive to the exact SFH and pro- 
duce populations of A^msp ~ 200 with an upper limit of 
A^MSP ^ 2000 for a range of reasonable SFHs. 
Upper limits on the number of CPs in the inner parsecs 
of the GC may be set using the detection of ^20 PWNe 
within 20 pc of Sgr A*. Using th e ~20 PWN candidates 
compiled by iMuno et al.l (|2008| ) in a catalog of diffuse 
X-ray sources and assuming PWNe are produced at a 
constant rate, we find that as many as 400 CPs may 
reside in the inner 20 pc. If the CP distribution follows 
that of the PWN candidates, then as many as 80 CPs 
could reside in the inner parsec. 

Finally, we consider measurements of the soft (kT w 
0.8 keV) component of the diffuse X-ray plasma in the 
inner 20 pc. Assuming that the plasma was heated by 
the transfer of kinetic energy from SNc, we may set 
an upper lirnit on the number of pulsars in this region 
(|Muno et al.ll2004l) . By considering two different cooling 
regimes, we find that the most extreme cooling scenario 
will produce a population of up to ^lO'^ CPs and '^10^/^ 
MSPs. Thus, we can set upper bounds of A^cp ^ 10'^ and 
A^MSP < 104(f^/10-2) for CPs and MSPs in the inner 
20 pc. 

Overall, we find that the above estimates are consistent 
with a CP population of Nqp < 100 and an MSP popu- 
lation of A^MSP ^ 10'^ in the inner parsec of the GC. 

10.2. Conclusions 

Current observations of the GC arc consistent with a 
population of up to ~'10"^ active pulsars beamed toward 
Earth within the central parsec around Sgr A*. This to- 
tal population may consist of up to A'^cp ^100 CPs and 
as many as A'msp ^ 10'^ MSPs. Such a population could 
distort the low- frequency measurements of both the spec- 
trum of Sgr A* and free-free absorption along the line of 
sight of Sgr A*. However, even with a potentially sizeable 
collection of pulsars, the difhcult observing conditions of 
the inner regions of the GC will make individual detec- 
tions a challenge. The strong interstellar scattering and 
large pulse broadening times mean that typical pulsar pe- 
riodicity searches at radio wavelengths will be sensitive 
to pulsars o nly at high o bserving frequencies and large 
bandwidths (jCordes fc Laziail997f) . Even if most of the 
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pulsars lie below the detection threshold, "giant" pulses 
intrinsic to the pulsar or as the result of an enhancement 
through multipath scattering may make a pulsar visi- 
ble a search for single pulses. Since there may be many 
more MSPs than CPs. one may search for pulsars us- 
ing interfcrometer imag ing surveys of compact radio ob- 
jects similar to that of iLazio fc CordesI (|2008[ ). Finally, 
search methods should be considered at wavelengths less 
affected by the scattering effects of the ISM. Since many 
MSPs produce gamma-ray emission, a directed search 
of the GC region with the Fermi LAT could potentially 
detect a pulsar in a blind periodicity search. 
Despite the difficulties in finding pulsars in the GC, 
the detection of a pulsar orbiting Sgr A* with an or- 
bital period of Porb ^ 100 yr would provide an un- 
paralleled test of gravity in the strong-field regime and 
could potentially allow the measuremen t of the spin 
and quadrupole moment of the MBH (iPfahl fc Loebl 
2n0llLaguna fc Wolszczanlll997l : 11^ fc KoneikinllTooa 
Liu et al.ll2012| ). Additionally, the detection of even one 
pulsar in the inner few parsecs would provide an excellent 
probe of the magneto-ionic material, the gravitational 
potential, and the star formation history in the vicinity 
of Sgr A*. In light of the significant scientific rewards 
and a potentially sizeable target population, we strongly 
recommend continued pulsar searches in the GC across 
a wide range of wavelengths. 
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